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C\| Abstract 

We consider new-physics (NP) contributions to the decay t —> bbc. We parameterize 
the NP couplings by an effective Lagrangian consisting of 10 Lorentz structures. We 
show that the presence of NP can be detected through the measurement of the partial 
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width. A partial identification of the NP can be achieved through the measurements of 
a forward-backward-like asymmetry, a top-quark-spin-dependent asymmetry, the partial 
rate asymmetry, and a triple-product asymmetry. These observables, which vanish in the 
standard model, can all take values in the 10-20% range in the presence of NP. Since 
■^j- |VtfoT4b| — iVts'Vcsl, most of our results also hold, with small changes, for t — > ssc. 
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I. INTRODUCTION 

On the whole, measurements of observables in the B system agree with the 
Standard Model (SM). However, some cracks have started to appear. There are 
now several quantities whose measured values differ from the predictions of the SM. 
Although these disagreements are not statistically significant - they are typically at 
the level of ~ 2c - they are intriguing since there are a number of different B decays 
and effects involved, and they all appear in b — > s transitions. Because of this, there 
have been numerous papers examining new-physics (NP) flavour-changing neutral- 
current (FCNC) contributions to the various b — > s processes. These analyses have 
been performed in the context of specific NP models, or model-independently 

In general, such NP can also contribute to FCNC processes involving the top 
quark. This has been looked at, though much less so than in B decays. However, 
given that the LHC will produce a large number of top quarks and will be able to 
measure flavour-changing t decays, it is important to explore the possibility of NP 
contributions to FCNCs in the top sector. In the past, analyses have focused on 
rare top decays such as t — > cV (V = g, 7, Z) and t — >• ch [HE]- Other top decays 
where NP effects have been examined include t — > br + u [2H9] and t — > W + dk [TO] . 

In this paper we examine the decay t — > bbc. In the SM, this decay occurs 
at tree level, via t — ¥ bW — > bbc. On the other hand, because it involves the 
small element V^ (~ 0.04) of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing 
matrix, the amplitude for this process is also rather small, and is therefore quite 
sensitive to NP. For example, there could be NP contributions to this decay in 
models with a charged Higgs boson [2] , or via t — > X°c — > bbc, where X° corresponds 
to some neutral particle (such as a Z' or a non-SM Higgs boson). Such processes 
could interfere with the SM process, leading to observable consequences, even if the 
intermediate NP particle were heavier than the top quark. Rather than restricting 
our attention to any one particular model, we examine NP contributions to t — > bbc 
model-independently (i.e., using an effective Lagrangian). 

Our model-independent treatment of t — > bbc takes into account the effects of the 
10 possible four-Fermi operators. These operators contribute to both CP-conserving 
and CP-violating observables. For the CP-even observables, we consider the CP- 
averaged partial width, a forward-backward-like asymmetry, and an asymmetry that 
depends on the spin of the top quark. For the CP-odd observables, we note that 
the decay t — ¥ bbc is dominated by one amplitude in the SM; i.e., there is only one 
weak phase involved. As such, all CP-violating asymmetries are very suppressed in 
the SM, so the observation of a non-zero asymmetry would be a smoking-gun signal 
of NP. In this paper, we discuss two types of CP-odd asymmetries: the partial-rate 
asymmetry (PRA) and a triple-product asymmetry (TPA). 

PRAs require a strong phase in order to be non-zero. Strong phases can arise 
due to gluon exchange, but it is expected that such phases will be small since the 
energies involved are so large. Another source of a strong phase is the width of 
the W. In our calculation, we ignore QCD-based strong phases and assume that 
the required strong phase is due entirely to the width of the W. This means that 
only SM-NP interference can lead to a PRA. On the other hand, in contrast with 
PRAs, TPAs do not require a strong phase in order to be non-zero. Thus, NP-NP 
interference terms can give rise to TPAs. As we will see, TPAs generated by SM-NP 



interference tend to be small, but NP-NP TPAs can be large. These are particularly 
interesting. 

We show that the measurement of the partial width by itself can reveal the 
presence of NP. However, if the NP exists, we will want to know its identity, i.e., 
which of the 10 operators is responsible, and the partial width measurement does 
not give us this information. In order to do this, it is necessary to measure the other 
quantities mentioned above. Since the various observables depend differently on the 
operators, the knowledge of their sizes will give us an idea of which NP operators are 
present. This will allow us in turn to deduce which model(s) might be responsible 
for the observed effects. 

Although we confine our attention to t — ¥ bbc in this work, we note that most 
of our results are easily transferable to t -> ssc by the replacement (b, b) — ¥ (s, s) 
in Feynman diagrams and expressions. Since |V«>\4&| — |Vt s V^ s |, the branching ratio 
for t — ¥ ssc is similar to that for t — ¥ bbc, and the two processes would apriori have 
similar sensitivities to NP effects. One difference between t — ¥ bbc and t — ¥ ssc is 
that the "CPT" correction to the PRA would be signific ant for t — ¥ ssc, whereas it 
is miniscule for t — ¥ bbc (see the discussion in Sec. IV A and Appendix [Bj). 



The remainder of this paper is organized as follows. In Sec. [TT] we write down 
the SM contribution to t — ¥ bbc, and also parameterize NP contributions to this 



decay in terms of an effective Lagrangian containing ten terms. In Sec. Ill (CP-even 
observables) we compute the CP-averaged partial width for the decay under consid- 
eration, as well as a forward-backward-like asymmetry and an asymmetry based on 
the spin of the top quark. The latter two asymmetries are both constructed in such 
a way that they are zero within the context of the SM. We close this section with 
a brief numerical study, noting that the CP-even asymmetries could reasonably be 



of order 10's of percent. Section [TV] contains our analysis of two CP-odd observ- 
ables - the partial rate asymmetry and the triple-product asymmetry. Section [V] 
concludes with a discussion of our results. Appendices [A] [B] and [C] contain some 
technical details. In particular, Appendix [B] contains a discussion of results related 
to the CPT theorem, namely which vertex-type corrections must be considered in 
the calculation of PRAs in order not to violate CPT. 



II. STANDARD MODEL AND NEW-PHYSICS CONTRIBUTIONS 

In this section we parameterize the NP contributions to t — ¥ bbc in terms of an ef- 
fective Lagrangian. We then write down expressions for the SM and NP amplitudes. 
These expressions are used in following sections to determine various CP-even and 
CP-odd observables. 

Figure hia) shows the Feynman diagram for the SM contribution to t — ¥ bbc. The 
resulting amplitude is given by, 

M w = -2V2G F m 2 w V cb V tb {u bla P L u t ) {u clp P L v b ) [-g a ^G T (q 2 )] , (1) 

where V is the CKM matrix. We work in the standard representation of the CKM 
matrix, in which V cb and Vt b are both real. Note that colour indices have been sup- 
pressed. The expression in square parentheses is proportional to the W propagator, 





(b) 



FIG. 1: Feynman diagrams for t —> bbc. Diagram (a) shows the SM contribution. Diagram 
(b) shows the NP contributions in the effective theory. The NP contributions are assumed 
to have the same colour structure as that of the SM. See Appendix [C] for comments 
regarding the more general case. 



with q = Pt - Pb = Pb + p c , G T {q 2 ) = [q 2 - rn 2 w + ie T {q 2 )}~ and e T {q 2 ) ^ q 2 T w /m w , 
where T w ~ 3G F m^/(2 v / 2vr). 

(Note: throughout this paper we neglect the leptons' and light quarks' masses. 
However, if this not done, the W propagator is modified to 
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G T (q 2 ) + q -4rG L {q 2 ) 
q 2 



(2) 



where Gl = \f^w + ^L^q 2 )} ■ ^r(q 2 ) and £h{.q 2 ) are related to the transverse and 
longitudinal widths of the W [6] , and they both depend on the quark masses. There 
has been considerable discussion in the literature regarding the correct form of the 
W propagator. (See, for example, Refs. [SHU ITH - H3] .) The above expression has 
been derived by performing a Dyson summation of the absorptive parts of the W 
self-energy diagrams in unitary gauge, with quarks and leptons in the loops. Some 
of the disagreement in the literature has focused on the form of Gl (see the brief 
discussion in Ref. [H], for example). Still, when all light masses are neglected, none 
of the observables in the present work depend on Gl- There seems to be broader 
agreement on the form for Gt in the literature, although many authors drop the q 2 
dependence in ex- Finally, we should note that the Pinch Technique may be used to 
reorganize perturbative calculations - even those involving resonances - in such a 
way that results are explicitly gauge-invariant (see, for example, Ref. [IS])- Rigorous 
application of the Pinch Technique to the problem at hand is beyond the scope of 
this work.) 



We parameterize new-physics effects via an effective Lagrangian C eS = C v 



£fff + £rfF» where ^ 



cff 



c * = j^{n v LL h,p L tcrPLb + n v LR h,PLtcrp R b 



+U V RL bjpPntcrfPLb + TZl R b^P R tcYP R b} + h.c, (3) 

£f ff = jf 2 {K S LL bP L t cP L b + U S LR bP L t cP R b 

+n s RL bP R tcP L b + n s RR bP R tcP R b} + h.c, (4) 

ATs = |tj {Crba^tca^b + iC TE ba^tca^be^^} + h.c. (5) 



In the above expressions, g' is assumed to be of order g, M is the NP mass scale and 
the TZ and C couplings may include weak (CP-violating) phases. For the Levi-Civita 
tensor, we adopt the convention e 0123 = +1. The NP contributions to t — ¥ bbc are 
illustrated in Fig. fltb). Colour indices are not shown in the above expressions, but 
are assumed to contract in the same manner as those of the SM (i.e., b with t and 
c with b). In FCNC models - those with a flavour-changing neutral particle such as 
a Z' or a scalar - the colour indices would contract in the opposite manner (i.e., c 
with t and b with b). It is straightforward to incorporate colour-mismatched terms 
into the effective Lagrangian. This topic is discussed further in Appendix [C} 
It is useful to define 

v _(9'\ 2 (rnwy n v LL = V2 g' 2 U V LL 

LL \g) \MJ V tb V cb 8G F M 2 V tb V cb } [) 

and similarly for the other 1Z and C couplings. In terms of the "X" parameters, we 
have the following expression for the NP contribution to t — > bbc, 

M NP = AV2G F V cb V tb {X V LL u bl ^P L u t u cl »P L v b + ...}. (7) 

To get a sense of the possible order of magnitude of the X couplings, note that, if 
g' ~ 2g and M ~ 500 GeV, then X\ L ~ 2.5 x K V LL . Thus, X V LL could reasonably 
be assumed to be of order unity. In other words, the SM and NP contributions 
to t — > bbc can very well be about the same size. When computing the effect of 
NP on a particular observable, it is therefore important to include both the SM-NP 
interference and NP 2 pieces. 

At present, there are no direct constraints on the X couplings. The precision 
measurements of V cb place an indirect constraint via the loop diagram shown in 
Fig. [2] It is known that some care must be taken when attempting to incorporate 
terms from an effective Lagrangian into loop calculations [16]. Incorporating the 
diagram shown in Fig. |2j we find the following expression for the effective Lagrangian 
for b — > civ, 

£ sm+np ^ _ 2 ^2G F v cb [(1 + Ql L ) c Ll ^b L + £l R CRlM h-fvL + h.c, (8) 




FIG. 2: Loop- level contribution of the NP operators to b — y civ. This contribution affects 
the measurement of V c b- 



in which we have dropped corrections of order O (mb/m t ). Using the Feynman rules 
for the various vertices, and employing dimensional arguments, we estimate 

Cll(R) ~ 9 /9 2 ^ Vtb ^ Xll{r) ' ^ 

Since semileptonic b — y c transitions are used to determine V c b, the experimental 
value of V c b can be used to bound X\ L and X\ R . Let us first consider the X\ L term 
in Eq. (J8J), ignoring the X\ R term. The X\ L term has exactly the same structure as 
the SM term. Its effect is thus simply to multiply any inclusive or exclusive b — y civ 

width by a factor of 1 + 2Re (Cll) + |Cll| • The current experimental value of V c b 

is V c b = (40.6 ± 1.3) x 10~ 3 \T7\, implying a 6.4% uncertainty on Vj. If we assume 
that the X\ L contribution to b — y civ is hiding in the experimental uncertainty of 
V^, we find the bound, 

Re (X V LL ) < 2.6, (10) 

in which we have neglected the quadratic contribution of (Yl> since it is small. 
Since XY R is associated with the right-handed quark current in Eq. QSJ) , its effect 
is process-dependent. For example, for B — y Div, the hadronic matrix element is 
only sensitive to the vector part of the hadronic current, so left-handed and right- 
handed couplings both have the same effect, and they can be absorbed in with the 
SM current [IS]. For other modes, such as B — y D*iv, the right-handed and left- 
handed currents must be treated differently |18| . Since our expression in Eq. (J9l) is 
somewhat of an approximation in any case, we assume the same upper bound for 
XY R as for X% L , i.e., 

Re (XY R ) < 2.6. (11) 



The NP operators considered in this work could contribute, via loops, to other 
observables as well. As an example, consider the decay B — > ipK s , which proceeds 
at tree-level in the SM. In the present context, the NP operators contribute to this 
decay via a diagram similar to Fig. [2| but with Vi replaced by cs in the final state. 
The resulting effective Lagrangian would be very similar to Eq. (J8J), which could 
lead to effects in the measurement of sin(2/3) [19]. We do not consider such effects 
further here. 

Although we do not perform any model calculations in this work, it is worthwhile 
to consider which types of models could give rise to the various NP operators. 
The terms in Eqs. (pi) and Q arise in models that contain new charged vector or 
scalar bosons. For example, extensions of the SM containing gauge bosons with left- 
and right-handed charged-current couplings (such as the Left-Right Model) would 
contribute terms such as those appearing in C\ - including the TZlr an d ^Vrl 
terms if there were some amount of mixing between the left- and right-handed gauge 
bosons. Models containing charged scalars (such as the charged Higgs bosons that 
appear in many extensions of the SM) could give rise to the terms in the expression 
for C^ s . Alternatively, there are many FCNC models containing a heavy neutral NP 
particle (such as a Z' or a neutral Higgs boson) with flavour-changing t-c couplings. 
Here Fierz rearrangements of the eight operator combinations {'J^Pi^r) [7 m -Pl,r] and 
(Pl,r) [Pl,r] (in the notation employed in Ref. [20]) lead to all ten of the operator 
combinations in Eqs. (p])-(|5]). In this case there would be mismatched colour indices 
between the NP and SM diagrams. Appendix [C] explains how to deal with this 
situation. 

In the following sections we compute various observables associated with the 
decays t — > bbc and t — > bbc. We take as our starting point the expressions for the 
SM and NP amplitudes in Eqs. ((!]) and ([7]), respectively. 

III. CP-EVEN OBSERVABLES 

We begin by considering three CP-even observables associated with the decays 
in question. The first is the CP-averaged partial width, normalized to the SM 
result; the second is a forward-backward-like asymmetry; and the third is a CP- 
even asymmetry that employs the spin of the top quark. In Sees. |IH AjjTH~C] we 



work out expressions for the various observables. Section III D contains a numerical 
analysis and discussion of the results. 

A. CP-averaged partial width 

We first consider the partial width for t — > bbc. The expression for the differential 
partial width , inc luding the various NP terms from the effective Lagrangian, may be 
found in Eq. ( Al ) in Appendix LAi Performing the integrations over p 2 = (p t — p c ) = 
(Pb + Pb) and q 2 , we find, 

- - f AY 

r (t -> bbc) ~ r SM (t -)■ bbc) I 1 + —^ [-0.04 x Re {Xl* L ) + Im (X]£)] 



+ 



+ 



3G F m 2 t 



y/2-K* (1 - <&)* (1 + 2Qv> v 



\YV I 2 _l I vv I 2 



X 



RL\ 



X 



RR\ 



Ar.r.l + A. 



■it 



L LiJ 



+ X 



RL 



X 



flfl 



+ 24|X T r + 96|X TS | 



'121 



where (jy = mw/tTit and 



r SM (t -► &6c) ~ ^^ (^) 2 (i - C 2 wT (i + 2C^) • 



24\/2 



7T 



(13) 



In calculating the expressions in Eqs. (1121) and (13), we have used the narrow width 
approximation, in which \Gt{q 2 )\ 2 is replaced by a 5-function in q 2 , appropriately 



normalized. [One exception is the term proportional to Re(X££) in Eq. (12), which 
was computed numerically] 



(Note that the term proportional to Im(X££) in Eq. (12), which is involved in the 



partial rate asymmetry discussed below, is not quite complete. In its present form, 
it would lead to a violation of CPT. To avoid running into problems with the CPT 
theorem, certain vertex-type corrections need to be included in the calculation when 
computing partial rate asymmetries. We discuss these extra terms in Appendix |BJ 



The partial width for the CP-conjugate decay may be obtained from Eq. (12) by 
complex conjugating all weak phases. 1 This has the effect of changing the sign of 
the Im(X^£) term, while leaving the other terms unchanged. Adding the widths for 
t — y bbc and t — y bbc, and dividing by twice the SM result yields 



K 



r + r 
2Psm 

1 + 0.0845 x 



O.OSxRe (Xlt) + \X V LL \ 2 + \X V LR \ 2 + \X V RL \ 2 + \X 



V |2 
RR\ 



+ \ (|Xf L | 2 + (Xf^l 2 + |X| L | 2 + \X S RR \ 2 ) + 24 |X T | 2 + 96 \X TE \ 



(14) 



in which we have inserted the known values for the various physical constants, and 
used the expression for T\y noted below Eq. ([!]). Note that, in practice, the term 
proportional to Re(X^£) is always small compared to the other terms. 



Above, we noted that the X's could well be 0(1). Thus, from Eq. (14), we see 
that the CP-averaged partial width could be used as a tool to search for NP. In 
particular, an experimental value for 1Z that is different from unity would give clear 
evidence for NP. (Depending on the size of the NP signal, it may be important to 



include corrections to Eq. (13) [21J.) 

On the other hand, all 10 NP operators contribute to 1Z in similar ways. Thus, 
even if the measurement of TZ reveals the presence of NP, it does not give us any 
clue as to the type of NP. For this reason, it is important to look for signs of NP 
in other quantities. This is most easily done using observables that are strictly 



1 Note that the "i" multiplying Xte in the NP amplitude docs not get complex conjugated when 
computing the amplitude for the CP-conjugate process. 
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zero within the context of the SM. Such observables would typically depend upon 
differing combinations of NP parameters, so that the observation of NP effects using 
several different observables would yield insight into the precise nature of the NP. 
In the following, we construct several asymmetries that are zero within the context 
of the SM and discuss their potential usefulness as tools for searching for NP. 

Note that the ratio TZ, defined above, will also appear in the denominators of all 
asymmetries considered below. Since TZ is primarily a sum of positive quantities, 
and since it will always appear in the denominators of the asymmetries, it will 
always tend to decrease the values of the asymmetries compared to the analogous 
expressions employing the approximation 1Z « 1. 



B. Forward-Backward-Like Asymmetry 

A tool that has historically been useful to experimentalists is the forward- 
backward (FB) asymmetry. The differential width for t -> bbc may be written 
in terms of q 2 = (p t — p^) = {pi + p c ) and cos 9, where 9 is the angle between the 
momentum of the top quark and that of the charm quark in the b-c rest frame. The 
FB asymmetry makes use of the following asymmetric integration over cos 9, 



FB 



dT 



-dcos 9 



dT 



-dcos 9 



dq 2 . 



i dq 2 dcos 9 J_ 1 dq 2 dcos 9 

We choose instead to work with the kinematical variables q 2 and p 2 , noting that 

^ 2 2p 2 - q 2 



(15) 



cos 9 



m: 



<f 



m. 



(16) 



Using Eq. (|16|), we may rewrite Eq. (15) as follows, 

dT 



FB 



mi-q" 



(m?-ff»)/2 d 1 2d P' 



-M 



" 2 V 2 dT 



dq 2 dp- 



-dp 2 



dq 2 



(17) 



Let us first consider the SM contribution to the FB asymmetry. The SM-only 
contribution to the width is such that 



dT 



SM 



dq 2 dp 2 



oc \G T (q 2 )\ 2 (q 2 + p 2 ) (m 2 - q 2 - p 2 ) 



see Eq. (Al) in Appendix Ah. Using the integration prescription in Eq. (17), and 



assuming the narrow-width approximation, we find, 



4SM 
^FB 



r SM , p SM 
1 FB t L FB 



rsM + r 



SM 



pSM 
1 FB 

Tsm 



3C 



w 



2 (1 + 2C 2 



-0.228. 



(19) 



W J 



Thus, we see that the SM contribution to the FB asymmetry is non-zero. 

As noted above, in order to use a particular asymmetry as a discriminator of NP, 
it is useful if the asymmetry is zero when no NP contribution is present. The regular 
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FIG. 3: Phase space for t — > bbc. The gray vertical bar shows the location of the W 



resonance at q 



«Ti2» v 



rriyy The shaded and clear regions (separated by the curve denoted 
"pr" ) are used in the construction of the FB-like asymmetry A-^2. The line indicated by 
"pp B " shows the boundary used in the usual definition of the FB asymmetry. 



FB asymmetry does not satisfy this requirement, as is evidenced by Eq. (19). It 
turns out, however, that if we modify the p 2 integration prescription somewhat, the 
SM contribution can be made to disappear. That is, instead of breaking up the 
integral over p 2 at the point p| B = (m 2 
boundary to the value p 2 , 



q 2 ) /2, as is done in Eq. (17), we move the 



r> 



dT 



dq 2 dp 



;dp 2 



dT 



dq 2 dp 



-dp 2 



dq 2 



in which p 2 is chosen such that [see Eq. (18)] 



(q 2 + p 2 ) {m 2 t -q 2 - p 2 ) dp 2 



2 2 



Then, by construction, 



^SM 



m '- q2 dT 



SA1 



dq 2 dp' 



dp 2 



dT 



SM 



dq 2 dp- 



dp 2 



dq 2 = 0. 



(20) 



(q 2 + p 2 ) (ml - q 2 - p 2 ) dp 2 . (21) 



(22) 



The new boundary, p 2 , is g 2 -dependent and can be solved for numerically. 2 Figure || 
shows the phase space available for t — > bbc and also indicates the two boundary 



2 The equation for p 2 is cubic and can also be solved analytically, although the resulting expressions 
for the roots of the equation are not particularly enlightening. 



10 



choices described above. The vertical band indicates the location of the W resonance 
for the SM contribution. 

Equation (Al ) in Appendix [A] gives the expression for dT/dq 2 dp 2 . The p 2 depen- 



dence of the SM piece is given by (q 2 + p 2 ) (m 2 — q 2 — p 2 ); the SM-NP cross terms 
and the |A^ L | and pfjyj terms have this same p 2 dependence. Since the inte- 
gration prescription described above is engineered to eliminate the SM term when 
integrating over p 2 , these latter terms also disappear upon integration over p 2 in 
this manner. Performing the integration numerically for the other terms yields the 
following, 



r> 



3G 2 F m 5 t (V tb V cb 



,2 r 



167T 3 



0.155 [\X V |2 ' ]vV '-' 



^LR\ 



+ \x 



RL\ 



'S |2 



+0.0208 (\XZ L \~ + \X s RR f + |Xf K | 2 + \X S RL \ 2 
+0.310 Re [X T (X s L l + X S RR ) - 2X TE (X S L * L - X S RR )] 

+ 1.8l(|X T | 2 + 4|X TB | 2 ) 



(23) 



The above expression is CP-even; i.e., the analogous expression for t 
same. Finally, we form an FB-like asymmetry as follows, 



bbc is the 



^ 



iy + r^ 

r + r 



i 



0.0393 [\X V |2 ' lvV r> 



-LR\ 



+ \x 



RL\ 



+0.00528 [\X?r\ + \X 



^RR\ 



X 



LR\ 



-l\Y s I 2 

+ \ a rl\ 



+0.0786 Re [X T (X S L * L + X S RR ) - 2X TE (X S L * L - - X S RR )} 

+0.460(|A T | 2 + 4|A TS | 2 ) 



(24) 



By construction, this asymmetry is only non-zero if NP contributions are present. 



Section HID contains a discussion of the range of sizes that are possible for the 



FB-like asymmetry. At this point we note only that asymmetries of order tens of 
percent are possible. 

C. CP-even Spin Asymmetry 



The final CP-even observable that we consider depends on the spin of the top 
quark [22] • We construct this asymmetry in such a way that it will be zero in the 
SM and potentially non-zero in the context of NP. Equation (A2) in Appendix [A| 
contains the expression for the absolute value squared of the total amplitude, keeping 
only those terms that contain the top-quark spin four-vector. The term proportional 
to \Gt\ in that expression is the SM contribution. The next term [proportional to 
KeyGTXYl)] arises from the interference of the SM contribution with one of the NP 
terms. The remaining terms are purely NP in origin. Inspection of Eq. ( |A2[ ) reveals 
that the SM term is proportional to p^ ■ s t . (This is related to the fact that, in the 
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SM, the spin of the top is in the direction of the momentum of the b in the top's 
rest frame [23]-) Working in the top rest frame, we define 

%± = ±— -E~ ("c - %cos% c ) , (25) 

sm u bc 

where n^ = Pb( c \/\Pb(c)\ an d where % c is the angle (assumed to be between and 
7r) between the three-momentum of the b and that of the c, in the top's rest frame. 
The cosine and sine of this angle are given, respectively, by, 

m 2 (p 2 - q 2 ) - p 2 (p 2 + q 2 ) 
C ° bc " {m\ - p 2 ) (p 2 + q 2 ) ' 

• Q 2m tv / 'p 2 q 2 (m 2 - q 2 - p 2 ) 

sm0 bc = 72 9\ / 9 i — V\ • ( 26 ) 

(m 2 - p 2 ) (p 2 + q 2 ) 

The vectors s\\± are in the decay plane and are perpendicular to p\ by construction. 
Setting s^± = (0, %±)> we then have p^ ■ St,± = 0. Thus, the SM contribution to 
the amplitude squared disappears for these orientations of the top quark's spin. We 
can thus construct an asymmetry based on this spin configuration that will be zero 
within the SM, making it potentially a sensitive probe for NP. We first define, 

rii^[r(% + )-r (%_)], (27) 

where the factor of "1/2" is to account for the average over the top quark's spins. 



Using Eqs. (A2), (25) and (26), and incorporating the integration over phase space, 
we obtain, 

r G 2 F mUVt b V cb ) , , .., .. Iv 
in = „ < A rD — A 



7Qn 2 . xr-LH -|-/,'L| 

• I\yS I 2 \vS I 2 , \vS I 2 \ Y S I 2 



+2 Re [X T (X[* L - X S RR ) - 2X TE (X S L * L + X S RR )] - 96 Re[X T X* TE ] j . (28 
Finally, summing over the process and the CP-conjugate process, we obtain, 

r ii+rn 



An 



r + r 

0.0607 f /. v |2 , vV ,2\ 1 



I 2 i __ q i 2 i __ q i 2 i _.<j 1 2 



, , I vV | z \Y V \ \ f I l^ s I I Y s I -I- I Y 5 I I T s i 

■^-1 ^| a m| ~ | A i?x| j ~ ^ 1J a ll| _ | A /y?| "+" | A £fl| _ I a rl| 
+2 Re [X T (XH - X S R * R ) - 2X TE (X[* L + X S RR )] - 96 Re[X T X* TE ] } , (29) 

where we have used the fact that 12 v / 2G F m 2 /[357r(l - (^) 2 {l + 2^)] ^ 0.0607. 
A discussion of numerical values obtainable for the CP-even single-spin asymmetry 
follows in the next subsection. 
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TABLE I: Some representative values for the FB-like asymmetry A-^. and the CP-even 
spin asymmetry An. The value for 1Z is also included for each case. 



yV yV 

^LR ^RL 


vS vS 
^LL ^LR 


^RL ^RR 
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Xte 
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A t 
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7.4% 
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2.5 2.5 






1.3 


5.2% 


15% 




1 


-1 


0.5 


0.125 


1.7 


6.9% 


-14% 




1 


-1 


0.5 


-0.125 


1.7 


12% 


29% 




-2.5 


2.5 




0.25 


1.8 


21% 


0% 








0.5 


-0.25 


2.0 


11% 


36% 


2.5 2.5 










2.1 


24% 


0% 








1 




3.0 


15% 


0% 




2.5 


2.5 


1 




3.3 


28% 


0% 



D. Discussion of CP-even Observables 



In this section we have described three observables that are even under CP and 
that could be used to detect the presence of NP in the decays t — >■ bbc and t — > 
bbc. Should NP be discovered, detailed analysis of such observables could allow 
experimentalists to map out the nature of the NP. 

The first observable considered in this section was a ratio, 7Z, which was defined 
to be proportional to the CP-averaged partial width. Of the observables considered 
in this work, TZ would probably be the simplest to measure experimentally. Decisive 
experimental deviation from the SM value TZ = 1 would be evidence for NP. 

On the other hand, TZ cannot be used to distinguish the different NP operators. 
To do this requires the use of other quantities. It is useful to employ observables that 
give a null signal within the context of the SM. For such observables, a significant 
departure from zero would be a "smoking-gun" signal for new physics. In addition, 
since they depend differently on the various NP operators, the observation of non- 
zero values for these observables would help in identifying the type of NP. 

The usual forward-backward asymmetry for t — > bbc is expected to be non-zero 
within the context of the SM. It is possible, however, to alter the kinematical weight- 
ing that is used in defining the FB asymmetry in such a way that the resulting "FB- 
like" asymmetry is zero within the context of the SM. Equation (24) defines the 



FB-like asymmetry A-^. in terms of an asymmetric integration over the kinematical 
variable p 2 . The integration is engineered in such a way that the SM contribution 



disappears kinematically. In Eq. (29) we formed a CP-even asymmetry using the 



spin of the top quark. This asymmetry was also defined in such a way that it was 
zero within the context of the SM. A non-zero experimental signal for either of these 
asymmetries would indicate the presence of NP. 

Table |T| contains some representative values for the CP-even asymmetries A-pi and 
A\\, along with the corresponding value for the ratio 1Z in each case. The entries 
in the table are ordered from smaller 1Z values in the top rows to larger ones in 
the bottom rows. As is evident from the table, when TZ is close to unity (meaning 
that it may not be a very clear discriminator of NP), it is still possible to have 
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FIG. 4: Scatter plots of the CP-even asymmetries A^. and Ah for various combinations of 
the NP parameters. 



CP-even asymmetries that are on the order of several percent. For larger TZ values, 
the asymmetries A-^i and A\\ can reach into the 10's of percent. Note, however, that 
there are some NP scenarios in which Ay suffers cancellations or is zero, even if 
the NP parameters themselves are non-zero. For example, if A^J = AJ£ L , the 
contributions from these two parameters cancel in A\\. 

The CP-even observables are displayed in another manner in Fig. [4j which shows 
scatter plots of Aj2 versus 1Z and A\\ versus 1Z. The points in this plot were obtained 
by generating real random numbers for eight of the ten NP parameters over various 
ranges. [X\ L and X\ R were excluded, since they do not contribute to the numerator 
of either asymmetry.) Asymmetries were only plotted if 72. < 3. Again, it is evident 
that CP-even asymmetries of order a few 10's of percent are possible. 



IV. CP-ODD OBSERVABLES 

In addition to the CP-even observables considered in the previous section, it is 
also possible to construct CP-odd observables related to the decay t — > bbc. In this 
section we consider two such observables. The first is the partial rate asymmetry, 
which compares the partial width for the process to that of the CP-conjugate process. 
The second is a triple-product asymmetry, which is formed using the spin of the 
top quark and the three-momenta of two of the final-state quarks. To be non- 
zero, both of these asymmetries require the presence of at least two contributing 
amplitudes with a non-trivial relative weak phase. Let us first consider the partial 
rate asymmetry. 



A. Partial Rate Asymmetry 

The SM amplitude for t — > bbc is dominated by a single contribution. As such, the 
partial rate asymmetry vanishes. In the presence of NP, the partial-rate asymmetry 
(PRA) can be nonzero if there is a NP contribution to the decay with a relative 
weak phase. As can be seen in Eq. (12), there is one important NP piece of this 



type - X^ L . The contribution to the PRA then comes from the interference of the 
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SM jy-exchange amplitude with the X\ L term. What we see in this subsection is 
that the PRA can actually be of order several percent if the Lorentz structure of 
theNP is (V-A)x (V - A). 

We have noted above that a non-zero PRA requires the interference of at least 
two amplitudes having a non-zero relative weak phase. Another requirement is that 
these amplitudes have a non-zero relative strong phase. Strong phases can come 
from the exchange of gluons, but they can also emerge from the imaginary parts of 
loop diagrams that do not involve gluons. In particular, if an exchanged particle in 
the process has a resonance, there is a strong phase associated with the width of 
that particle. Strong phases originating from particles' widths have been used to 
generate PRAs in many different systems, including t — > bbc [21 [21], t — >■ br + v [21 E]- 
El [H|, and various supersymmetric decays [2"5T - [2"T] . In the present calculation, the 
width of the W provides the required strong phase. This means that the PRA can 
only arise from SM-NP interference, since NP-NP interference terms do not have a 
relative strong phase. 



Using the expression in Eq. (12), and recalling that the analogous expression for 
the CP-conjugate process involves the complex conjugation of the weak phases, we 
immediately find the following expression for the partial rate asymmetry, 

^^4^M^)^x M ^). (30) 

As was noted above, the PRA requires the existence of a non-zero relative strong 
phase between interfering amplitudes. In this example, the strong phase is provided 
by the width of the W, which is the reason that the PRA is proportional to T^- 



Examination of Eq. (30 ) reveals that the best-case scenario for the PRA occurs when 
X\ L is purely imaginary, or nearly so, and all other NP coefficients are zero. In this 
case, 1Z ~ 1 + 0.0845 \XYl\ , and we find that the PRA is maximized when | XYl\ — 
3.44. The maximum possible PRA, obtained in this manner, is approximately 18%. 
As is well-known [23 [2H], the CPT theorem requires that we actually be a bit 
more careful when computing PRAs. In particular, invariance under CPT requires 
that the total width of the top be equal to that of the anti-top. Our result in 



Eq. (30) shows that, under certain circumstances, the partial width for t — )■ bbc is 
not equal to the partial width for t — > bbc. This necessarily implies that there must 
be compensating partial rate asymmetries in other top/anti-top decay modes such 
that the total top width is still equal to the total anti-top width. In order to respect 
CPT in this way, it turns out that we need to include another class of diagrams, 
shown in Fig. [5} These diagrams contribute to various top decay modes, inducing 
partial rate asymmetries in these modes in such a way that the total top width is 
equal to the total anti-top width. In the special case t — > bbc, the effect is such that 
"IV" in the numerator of Eq. Q gets replaced by U T W -T (W ^bc)" [23 EH ED], 
which is to say that the strong phase due to the rescattering process W —)-bc —> W 
does not contribute to the PRA. Since T (W —¥ be) is very small, we may safely 
neglect its effect. It is worthwhile to explore this point a bit further, however, and 
we do so in Appendix [Bj Specifically, we verify that the diagrams in Fig. [5] interfere 
with their SM counterparts in such a way that the CPT theorem is respected, and 
we also comment on the PRAs that result in other decay modes due to the NP 
effective operators for t — > bbc. 
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FIG. 5: Vertex correction-type diagrams involving the effective operators shown in Fig. [T] 
(b). These diagrams contribute to the cancellations required by the CPT theorem. The 
dashed line indicates that only the absorptive parts of the diagrams are computed. 



B. Triple Product Asymmetry 

Mathematically, triple-product asymmetries (TPAs) in t — >■ bbc are related to 
terms of the form Vi ■ (vj x v k ) in the absolute value squared of the amplitude, where 
each of the Vi could represent a momentum or spin. Working in the rest frame of the 
top quark, there are only two independent three-momenta. Thus, in order to obtain 
a non-zero TPA, we need to include one or more spins in Vi ■ (vj x Vk). Since the light 
final-state quarks hadronize, it is difficult to gain useful information from their spins. 
The situation is different for the top quark, however, since it decays too quickly to 
hadronize. In this case, we can construct asymmetries based on s- (pi x p 2 ), where s 
is the top quark's spin and p% and p 2 are two of the final-state momenta |32j . In the 
context of the calculation, these terms arise from expressions such as e Q ^p"sfp^pf. 

Now, the PRA considered above contained a factor of Y\y ~ 2 GeV in the numer- 
ator. The presence of this factor was due to the requirement that there be a relative 
strong phase between diagrams contributing to the PRA. On the other hand, TPAs 
do not require a strong phase and are thus not suppressed by a factor of Yw This 
means that TPAs could in principle be much larger than the PRA considered above. 
As we shall see, there are in fact certain NP operators that can produce a large TPA. 

Because TPAs are CP-odd quantities, they require a non-zero relative weak phase 
between interfering diagrams, just as the PRA did. But because no strong phase is 
necessary, TPAs can in principle arise both from SM-NP and NP-NP interference. 
(Due to the strong phase requirement, the PRA could only arise from SM-NP in- 
terference.) What we find, however, is that the only TPA that survives is one due 
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TABLE II: Some representative values for the triple-product asymmetry. The second to 
last column also shows 1Z for each case. 
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X RR 


^Ty 


Xte 
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4 TP 
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-1.5a 


0.5 





1.6 


23% 
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1 
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to NP-NP interference. 

All triple-product terms in the absolute value squared of the amplitude may be 
written in terms of e a ^5Pt s tPiPc- Keeping only such terms, we find the following 
expression in the rest frame of the top quark, 

li: e witp 

colours b,b,c spins 

= l536G F m 2 t {V tb V cb f s- (p- b x p c ) Im[X T (X s L i+X s RR ) -2X TE (X S L * L -X S RR )} , (31) 



in which s denotes the top's spin [see also Eq. (A2)]. In computing the above 
expression, we have summed over quark colours and over the final-state quarks' 
spins, and have divided by 3 for the average over the top quark's colours. Setting 



we define 



s ±,± 



1 



Pb X Pc 

\Pb X Pel 



r TP = - [r ( S "i, + ) - r ( s v)] 



(32) 



(33) 



where the factor of "1/2" is to account for the average over the top quark's spins. 



Using the result in Eq. (31) and incorporating the integration over phase space, we 
obtain, 



II 2G F m t (VtbVeb) T ^ r v (yS*j_yS*\ ov ( vS* vS*X\ 
Itp - ^—^ lm [X T [X LL + X RR ) - ZX TE [X LL - X RR )\, 



357T 2 

in which we have used the fact that 

^ X & = i 
Finally, we define the TPA as 



r 2 2 1 i 



p 2 )] 



1/2 



4 TP _ 



r T p — r 



TP 



r + r 



so that 



/i CP _ 



1 A8V2G F m 2 t lm[X T (Xf£ + X S RR ) - 2X TE (X S L * L - X S RR )] 



K 35tt 
0.243 



;i-c^-) 2 (i + 2c 2 



WJ 



K 



lm[X T {XH + X S R * R ) - 2X TE {XH - X S RR )} . 



(34) 



(35) 



(36) 



(37) 
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FIG. 6: Scatter plot of the CP-odd TP asymmetry A^p f° r var i° us combinations of the 
NP parameters. 



Table [XT] contains some numerical results following from the above expression, show- 
ing that the TPAs can indeed be large - of order 10's of percent - if the NP coeffi- 
cients are assumed not to be suppressed. Figure [6] shows a scatter plot of A^p versus 
TZ. The points in this plot were obtained by generating combinations of purely real 
and purely imaginary random numbers for Xf L , X RR , Xt and Xte over various 
ranges. Once again, asymmetries were only plotted if 1Z < 3. It is evident from the 
plot that relatively large TPAs are possible. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we consider new-physics (NP) contributions to the decay t — > bbc. 
In the Standard Model (SM), this is a tree- level process: t — > bW — > bbc. However, 
the SM amplitude involves the small Cabibbo-Kobayashi-Maskawa element V c b (~ 
0.04), and is therefore suppressed. As a result, the decay is quite sensitive to NP. 
Rather than working within the context of any one particular extension of the SM, 
we parameterize the NP couplings by an effective Lagrangian that includes the 
10 possible four-Fermi operators. We show that the SM and NP contributions to 
t —7- bbc can indeed be about the same size. 

We first compute the t — > bbc partial width in the presence of NP. The ratio TZ, 



defined in Eq. (|14|), provides a quantitative measure of the deviation of the partial 

in the SM). This shows clearly that this 
significant deviations of TZ 



width from its SM expectation (TZ = 1 

observable is excellent for showing that NP is present 

from 1 are possible. 

On the other hand, the partial width is not a good observable to identify the 
new physics - all 10 NP operators contribute to TZ in a similar way. In order to 
get an idea of the type of NP present, it is necessary to consider other quantities. 
To this end, we construct two CP-conserving and two CP-violating observables: (i) 
CP-even: a forward-backward- like asymmetry \A-pi - Eq. (24)] and a top-quark- 



spin-dependent asymmetry [Ai\ - Eq. (29)], (ii) CP-odd: the partial rate asymmetry 
[^cp ~ Eq. (30)] and a triple-product asymmetry [A^p - Eq. (37)]. In each case, the 
observable is formulated in such a way that it is zero within the context of the SM. 
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The key point is that these observables depend on differing combinations of the NP 
parameters. This gives them different sensitivities to the various Lorentz structures 
present in the NP effective Lagrangian. 

The allowed values of these four observables vary greatly depending on the values 
of the NP parameters, but results in the 10-20% range are possible (see Tables ffl and 
[lT| ). If NP is present, it may well produce measurable values of these observables. 
Taken together, the measurements of these quantities will give a good indication of 
the type of NP present. 
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Appendix A: Useful Expressions for t — > bbc 



This Appendix contains two expressions that are used to compute observables 
in the main body of the paper. We take m^ = m c ~ 0. The first of these is the 
expression for the partial differential decay width for t — > bbc. Using Eqs. (hj) and 
(J7|) and averaging over the top quark's spins and colours, we find the following, 



dT 3G 2 F (V tb V cb ) 



dq 2 dp 2 ~ 2(2 7 r) 3 m i 3 



(g 2 + p 2 ) (m 2 - g 2 - p 2 ) 



mw 



V |2 
RR\ 



G T \ 2 + Am 2 w Re {G T X v L l) + 4 (| X V LL \ 2 + \X 

+ 4p 2 (m 2 -p 2 )(|AL| 2 +|X^| 2 ) 

+ g 2 K - g 2 ) (|X£ L | 2 + (A^l 2 + \X[ R \ 2 + \X 

+ 8 g 2 {-m 2 t + g 2 + 2p 2 ) Re [X T (X S L * L + X S R * R ) - 2X TE (X_ 



-S I 2 



s* 

LL 



X S nn)} 



+ 32 



m 



2 ( „2 



(g 2 + 4p 2 )-(g 2 + 2p 2 )^ (|A T | 2 + 4 \X TE \ 2 ) 



(Al) 



in which V t b and V c \, have been taken to be real. The analogous expression for 
t -¥ bbc is obtained by complex conjugating all of the NP coefficients [X\ L , etc.), 
while leaving Gt unchanged. 

It is also useful to have the expression for the absolute value squared of the 
amplitude, keeping only the terms that contain the spin four vector for the top 
quark. This expression is used to compute the CP-even single-spin asymmetry in 
Sec. Ill C and the TP asymmetry in Sec. IV B| Keeping only terms containing the 
spin four vector of the top quark, we find, 



\E E M 



colours bb c 



spins 
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9QG 2 F m t (V tb V cb f\ {m 2 t -q 2 -p 2 ) 



m 4 w \G T \ 2 -4m 2 w Re(G T XYt) 



\(\xU 2 -\xU 2 



P- b -St-4p 2 (\xV R \ 2 -\xV L \ 2 ) Pc .s t 

+ 8 Re [X T (Xf* - X* ) - 2X TE (Xf* + X*) ] [(m? - g 2 ) p s ■ « t + p 2 p b ■ S J 
+ 128 Re [X T X* TE ] [(2m 2 - q 2 - 2p 2 ) p> • s t + (q 2 + 2p 2 ) p c ■ s t ] 

- 16 Im [X T (X S L * L + X S R * R ) - 2X TE (X[* L - X S RR )] e {p t) s t ,p- h , Vc ) V (A2) 

in which s t denotes the top's spin four vector and e(p t ,s t ,p b ,p c ) = ^a/3-ysPt s tPj.P 5 c- 
In writing the above expression, we have used the fact that p t ■ s t = 0. We have also 
summed over quark colours and over the final-state quarks' spins, and have divided 
by 3 for the average over the top quark's colours. 

Appendix B: CPT and Beyond 

The CPT theorem requires the total decay width for the top to be equal to that 
for the anti-top. An apparent violation of the CPT theorem arises, however, if the 
NP contributions in Fig. [I] (b) are the only ones that are kept. That this is the 
case is straightforward to see, since the diagram in Fig. [I] (b) affects the partial 
widths for t — > bbc and t —¥ bbc differently [leading to the PRA in Eq. (30)], but 
has no effect on the other top or anti-top decay modes. Thus, the top and anti- 
top total widths are not equal if only such contributions are kept, resulting in an 
apparent violation of the CPT theorem. This phenomenon is well-known (see, for 
example, Refs. [T1I2" HI2"8"H3"T] ). In this Appendix we show that the inclusion of certain 
vertex-type corrections gives rise to compensating differences in the top and anti-top 
widths. The sum of the differences is zero, so that the top and anti-top widths no 
longer differ, in agreement with the CPT theorem. 

Let us define the partial width difference for the decay t — > bjk as follows, 

Ar (t -> b]k) = T (t -> b]k) - T (t -» bjk) , (Bl) 

in which j and k could refer either to quarks or to leptons. For the case t — > bbc, 
the main contribution to Ar is due to the interference between the SM and NP 
diagrams indicated in Fig. [Tl Another important set of contributions for the decay 
t — > bjk is indicated in Fig7[5J The absorptive parts of these vertex-like corrections 
interfere with their associated SM diagrams in such a way that the conservation of 
CPT is manifest. Using the Cutkosky rules to calculate the absorptive part of the 
vertex-like corrections, we find 

Ar (t -»- bjk) ~ 2V ^ GFrr ^ iVcb)2 lm(Xll)T (t -> bW) [5- fb 5 kc -B (W ->jk)}. (B2) 
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Summing over j and k (including both quark and lepton final states), we have 

J2&r(t^b]k)=0, (B3) 

demonstrating that the CPT theorem is indeed respected once the absorptive parts 
of the diagra ms i n Fig. [5] are included. 



Equation (B2) gives a correction to the PRA for t — ¥ bbc, leading to the following 



modification of Eq. (30) 



Im (X^*\ 
Acp~ 0.102 x \-^- [1 - B (W ->• 6c)] . (B4) 

The correction to the original expression is miniscule, since B (W — > be) ~ 
\V cb \ 2 /3 ~ 5.5 x 10" 4 . 

An interesting consequence of the CPT theorem is that, if NP operators give 
a PRA in a particular decay mode (such as t — > bbc, as in our case), then those 
same NP operators must also contribute to one or more other decay modes in such 
a way that the total width of the top is the same as that of the anti-top. This 
means that those other decay modes must also have PRAs (barring other accidental 



cancellations). We can use Eq. (B2) to estimate the PRAs in other decay modes 
due to the NP operators in Eqs. (J3j)- (Jsj) . The resulting expression is given by, 

A CP (t -> b]k) ~- ^ GFm2w{Vcb) \ m(xYl) --5.6 x 10- 5 Im(A^) , jA^6c.(B5) 

Thus the contributions of these NP operators to PRAs in other top decay modes 
are expected to be very small. One could also consider the complementary question: 
Are there NP operators, other than those given in Eqs. p])-©, that could contribute 
to the PRA in t — > bbc! The answer to this question appears to be yes. For example, 
the effective operators (sOic)(c02b) or {dOiu)(c02b) could appear in a diagram 
similar to that in Fig. [5j but with the usual SM tbW vertex at the top, and the 
NP-induced one-loop correction to the Wbc vertex at the bottom. Such operators 
are constrained by B decays. 

We should note that the NP effective operators in Eqs. ([3])-((5]) give rise to other 
loop-level diagrams that could contribute to PRAs in top decays. The contributions 
in different decay modes would still complement each other in the sense that the 
total top and anti-top widths would remain equal. Figure [7] shows an example 
of loop-level corrections to t — > ffc (/ = d, s, b) mediated by the NP operators 
considered in this work. These diagrams could interfere with their corresponding 
SM diagrams to induce PRAs. We do not compute such contributions here. 



Appendix C: Effect of Including Colour- mismatched Terms 

The effective Lagrangian incorporating NP effects given in Eqs. ([3])-((5| assumed 
that the colour indices contracted in the same manner as those of the SM diagram. 
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FIG. 7: A loop-level contribution of the NP operators that could contribute to PRAs in 
t — > ffc, with f = d,s, b. 
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This need not be the case, so it is useful to consider the effects of including colour- 
mismatched terms in the effective Lagrangian. To this end, let us generalize the NP 
effective Lagrangian in Eq. (pi) as follows, 

,12 

{n v LL b alyb P L t a c b YP L b b + K V L ' L b alfl P L t b c b YP L b a 

+K V LR b alii P L t a c b ^P R b b + U V L ' R b alil P L t b c b ^P R b a 

+ ...} + h.c, (CI) 

and similarly for Eqs. Q and (J5J). In this expression, the subscripts a and b are colour 
indices and the primed coefficients correspond to the new, colour-mismatched terms. 
The total amplitude for t a — y b b b c Cd (with the subscripts a, b, c and d representing 
the colours) could then be parameterized as 



M abcd = Y^ i^abScd + K'AAc) Mi 



(C2) 



in which the sum runs over the SM diagram, plus all NP contributions. The factors 
IRj and M/ t are the coefficients for the colour-matched and colour-mismatched terms, 
respectively, and are assumed to contain all of the weak phases. (The W coefficient 
for the SM term is assumed to be zero.) For a given value of i, the phases of M; and 
W t could be different. The factors Aii contain all the spinors and 7 matrices and, 
in the case of the SM diagram, the W propagator. 

Summing over the quarks' colours and dividing by 3 for the average over the top 
quark's colours, we find, 
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Im 



+ 
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+ (* 



+ 



r* 



* j 



i[R*) 



Re(MiM*) 



lm(MiM*) . (C3) 
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The IRjIR* terms in the above expression correspond to the "colour-matched" terms 
that we have taken into account in this work. The other terms are new. 



Equation (C3) can be used to generalize the expressions in this paper, provided 



the expressions have already been split cleanly into pieces containing the weak phases 
(Rj, etc.) and those containing the spinors and any strong phases (Mi). Expressions 
containing SM-NP cross-terms may safely set Msm = 1 and incorporate the entire 
amplitude into the "A4sm" part [in Eq. (C2)], since V t b and V c b have been taken to 
be real. As an example, the generalized form for Eq. (14) would be 



K ~ 1 + 0.0845 x 



0.05 x Re [XY? + -X 



Vi* 



+ \ a ll\ 



+ \X V L ' L \ 2 + \Re(Xl L Xl'*) + . 



(C4) 



in which we have used the fact that 
0.243 



°SM 



0. Similarly, Eq. (37) would become 
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CI' 



K 



Im{X T (X[ I + X S RR ) - 2X TE (X[* L 

I V' I vSl* i vSl*\ o Y"' ( vSi* 
+ A T l A LL + ^Rr) - ZJi TE (^LL " 

fy ( vSl* i vS>*\ O V (~V~S/* 

- [A T [Ji LL + JL RR ) - IJ^TE \^_ 



vS* 



X S rr) 



+ 



LL 



ItRJ 
RRJ 



ySf: 



vi I vS* i vS* \ o -yf I vS* vSt \11 

-a t [Ji LL + Ji RR ) - 2A TE [ji LL - yi RR )\ ) 



and Eq. (30) would become, 



.4 



CP 



0.102 

n 



x Im 



LL + 3 A LL 



(C5) 



(C6) 



Finally, an expression such as "^(GtXYI)" in Eq. (Al), which contains both a 
strong phase (in Gt) and a weak phase (in X\ L ), first needs to be separated into 
two pieces using Re (AB) = Re(A)Re(B) - lm(A)lm(B). 
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